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Pion-induced Drell-Yan processes and the flavor-dependent EMC effect 
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Pion-induced Drell-Yan processes are proposed as a potential tool to measure the flavor depen- 
dence of the EMC effect, that is, the flavor-dependent modification of quark distributions in the 
nuclear medium. Existing pionic Drell-Yan data are compared with calculations using a recent 
model for nuclear quark distributions that incorporates flavor-dependent nuclear effects. While no 
firm conclusions can yet be drawn, we demonstrate that existing Drell-Yan data seem to imply a fla- 
vor dependence of the EMC effect. We highlight how pion-induced Drell-Yan experiments on nuclear 
targets can access important new aspects of the EMC effect, not probed in deep inelastic scatter- 
ing, and will therefore provide very stringent constrains for models of nuclear quark distributions. 
Predictions for possible future pion-induced Drell-Yan experiments are also presented. 
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The depiction of the nuclear structure functions in the 
valence quark region was discovered in 1983 by the Eu- 
ropean Muon Collaboration (EMC) in a muon-induced 
deep inelastic scattering (DIS) experiment [l| , and is now 
known as the EMC effect. The EMC effect provides clear 
evidence that the quark distributions in nuclei are mod- 
ified compared to those of free nucleons. The observa- 
tion of the EMC effect was unexpected, since the nuclear 
binding energy is orders of magnitude smaller than the 
energy scale probed in DIS. Other DIS experiments using 
electron (2] , muon [l|, y[ and neutrino [4[ beams have now 
confirmed the EMC effect over a broad range of nuclear 
masses (A) and momentum transfers (Q 2 ). The improved 
precision of these experiments, including the recent ex- 
periments on light nuclei (a, lf| , have provided a detailed 
and multi-dimensional view of the nuclear modification 
of the quark distributions. The EMC effect has also been 
experimentally verified in the time-like region usingboth 
the pion- and proton-induced Drell-Yan processes [iHSl]- 

Despite a quarter century of significant experimental 
and theoretical effort, the specific origins of the observed 
A dependence of the nuclear quark distributions have yet 
to be unambiguously identified. Attempts to explain the 
EMC effect have led to a large collection of theoretical 
models [10, LLU] , many of which are capable of describing 
the essential features of the data, however the underlying 
physics mechanisms in each model are often very differ- 
ent. Therefore, it appears likely that before the true ori- 
gins of the EMC effect are understood new experiments 
are necessary, which explore aspects of the EMC effect 
not effectively probed in DIS, and will therefore help dis- 
tinguish between the various models. In this Letter we 
explore the effectiveness of pion-induced Drell-Yan pro- 
cesses as a tool to investigate new aspects of the EMC 
effect PJl, LL3( j m particular its flavor dependence. 

A new calculation of the modifications of nucleon quark 
distributions in the nuclear medium has recently been re- 
ported [lj, LL5[. In this approach, by Cloct, Bcntz and 
Thomas (CBT), the Nambu-Jona-Lasinio model is used 



to describe the coupling of the quarks in the bound nucle- 
ons to the scalar and vector mean fields inside a nucleus. 
These nucleon quark distributions are then convoluted 
with a nucleon momentum distribution in the nucleus to 
generate the nuclear quark distributions [14j . An impor- 
tant feature of this model is that for TV ^ Z nuclei (where 
N and Z refer to the number of neutrons and protons) 
the isovector-vector mean field (usually denoted by p°) 
will affect the up quarks differently from the down quarks 
in the bound nucleons. Therefore, this model has a novel 
prediction that the u and d quarks have distinct nuclear 
modifications for N / Z nuclei. In other words, the 
EMC effect is flavor sensitive and depends on the N/Z 
ratio of the nucleus |1J . This flavor-dependent EMC 
effect was also shown [15|, [l6[ to account for a large frac- 
tion of the dis crep ancy between the NuTeV collaboration 
measurement [171 ] of the weak mixing angle and the world 
average. 

Figure [T] shows the prediction for the nuclear depen- 
dence of quark distributions in the CBT model. The solid 
curve is the result for the usual EMC effect in symmetric 
(N = Z) nuclear matter, that is, the ratio FA/FJP , where 
FA and F/P are the per-nucleon structure functions of the 
nucleus and the deuteron, respectively. 1 The calculation 
is in good agreement with the data extracted by Sick and 
Day [18[ . For an isoscalar nucleus the u and d quark dis- 
tributions are modified in the same manner, however for 
& N =/= Z nucleus the u and d quark distributions are 
predicted to undergo different modifications in the CBT 
model. For example, the dashed and dot-dashed curves 
in Fig. [1] are, respectively, the predicted ratio of the u- 
and d-quark distributions in a gold nucleus to those in 



The CBT results used throughout this paper arc calculations 
performed for asymmetric nuclear matter. Results for particu- 
lar nuclei are obtained by choosing the N/Z ratio equal to the 
nucleus in question. The deuteron results are obtained from a 
combination of free proton and neutron quark distributions. 




FIG. 1: Ratios of quark distributions and structure functions 
in nuclear matter versus the deuteron plotted as a function 
of Bjorken-s, at Q 2 = lOGeV 2 . The solid circles are data for 
N — Z nuclear matter from Ref. [Ig|- The solid curve is the 
calculation of F 2 jF® for N — Z nuclear matter from Cloet, 
Bentz, and Thomas [El, [HI]. The dashed and dot-dashed 
curves are the ratios of quark distributions in a gold nucleus 
to those in a deuteron, for u and d quarks, respectively. 



the deuteron. For an N > Z nucleus (such as gold), the 
p° mean field leads to a stronger nuclear binding for the 
u quarks compared to that for the d quarks pjj]. Hence, 
the medium modification of the w-quark distribution is 
enhanced, while the d-quark distribution is modified less 
by the medium, as illustrated in Fig. Q] 

The flavor dependence of the EMC effect is a promis- 
ing experimental observable to distinguish among the 
plethora of models that can describe the EMC effect. 
In the simplest picture, we can consider the nuclear par- 
ton distribution functions (PDFs) to be just the nucleon 
PDFs smeared by the Fermi motion of the nucleons. In 
this picture it is natural that the nuclear modification 
would be very similar for u and d quarks. Similarly, 
the Q 2 -rescaling model [19| also predicts no quark-flavor 
dependence for nuclear modifications. In contrast, the 
pion-excess model j20l - |22l | would naturally predict flavor- 
dependent nuclear modifications arising from the differ- 
ent isospin composition of the pion cloud components of 
protons and neutrons. 

Since the inclusive DIS experiments measure the com- 
bined nuclear modifications of u and d quarks, it is diffi- 
cult to extract the quark-flavor dependence of the EMC 
effect in these experiments. However, it has been demon- 
strated [23[ that the parity violating DIS asymmetry, 
Apy, is very sensitive to the quark-flavor dependence of 
the EMC effect. Measuring Apy for a variety of N ^ Z 
nuclei has the potential to provide a clean demonstra- 
tion of flavor-dependent nuclear modification of the the 
quark distributions. Such a measurement is part of the 
large parity violating DIS program proposed [24j at the 
upgraded 12 GeV JLab facility currently under construc- 
tion. 

Semi-inclusive DIS (SIDIS) on heavy nuclear targets, 
in which the flavor of the struck quark is tagged by the 
detected hadron, is also a promising experimental tool to 



search for the flavor-dependent EMC effect. Recently, Lu 
and Ma [25[ pointed out that charged lepton SIDIS off 
nuclear targets and the deuteron can be used to probe the 
flavor content of the nuclear quark sea, which can help 
distinguish between the various models of the EMC ef- 
fect. Indeed, a SIDIS experiment [26| aiming at a precise 
determination of flavor dependence of the EMC effect has 
also been proposed at the upgraded 12 GeV JLab facility. 

The focus of this Letter is pion-induced Drcll-Yan pro- 
cesses, which are complementary to the DIS processes 
mentioned above and provide another experimental tool 
with which search for flavor-dependent effects in the nu- 
clear modification of the nucleon structure functions. 
First, we will identify the experimental observables sen- 
sitive to the flavor-dependent EMC effect. Then we com- 
pare existing pion-induced Drcll-Yan data with calcula- 
tions using the CBT model, and examine how well the 
data can constrain the flavor dependence of the EMC ef- 
fect. Finally, predictions for possible future pion-induced 
Drell-Yan experiments are presented. 

The leading order Drell-Yan cross-section for a pion 
interacting with a nucleus is given by 



dcr 2 ± , 



4na 2 



dx-n dX2 9 S X-„ X2 

X Z^2 6 9 feWivZ* ) Va(x 2 ) + q-K±(Xn) Qa(X2)] , (1) 
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where a is the fine structure constant, s is the center 
of mass energy squared, x v is the Bjorkcn scaling vari- 
able for the interacting quark in the pion and x 2 is the 
analogous quantity for the nucleon in the target nucleus. 
Quark flavor is labeled by q, where e q is the quark charge, 
the pion quark or antiquark distributions are labeled by 
q^± and q„± , respectively, and the subscript A indicates 
nuclear quark or antiquark distributions. 

To explore the sensitivity of pion-induced Drell-Yan 
processes to a flavor-dependent EMC effect, we consider 



the three ratios 
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where A represents a nuclear, D a deuteron and H a 
hydrogen target. Assuming isospin symmetry, which im- 
plies u OT + = d^- , u„- = d n +, u n + = (In- , u^- = d^+ and 
keeping only the dominant terms in each cross-section, 
one readily obtains 
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The up and down nuclear quark distributions are labeled 
by ua and d A respectively, ud is the up quark distribu- 
tion in the deuteron and u p the up quark distribution in 
the proton. Eqs. ©-(HJ) demonstrate that these Drell- 
Yan cross-section ratios are very sensitive to the flavor 
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FIG. 2: The existing data for the ratios 
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sus the predictions using nuclear PDFs of the CBT model for 
tungsten (red solid) and N = Z nuclear matter (blue dashed) . 



dependence of the EMC effect. Moreover, these ratios 
are not sensitive to the uncertainty of the pion structure 
functions, which are not yet determined accurately. 

To study the sensitivity of the pion-induced Drell-Yan 
processes to the flavor dependence of the EMC effect we 
calculate these ratios using the nuclear and nucleon PDFs 
from the CBT model [1J, [l5[ • Rather than using the ap- 
proximate expressions of Eqs. ©-((I]), the expression for 
the Drell-Yan cross-section in Eq. ([T]) is used in the cal- 
culations. Terms involving heavy quarks (s, c, b, t) are 
not included due to their negligible contributions. The 
pion-induced Drell-Yan cross-section for each nuclear tar- 
get is determined using the nuclear PDFs from the CBT 
model. The free nucleon PDFs obtained from the CBT 
model are used to calculate the Drell-Yan cross-section 
from the deuteron, with ujj = (u p + d p )/2 = do- For the 
pion PDFs, we use the GRV-P LO parametrization [271 ]. 
As expected, the results are found to be insensitive to 
the choice of the pion PDFs. 

In Fig. [2] we compare our calculations of the pion- 
induced Drell-Yan cross-section ratios with the existing 

data. The top left panel shows the ratio of ^dv ^-+pi 
from the NA10 experiment pj. These plots contain both 
the Pbeam = 286 and 140 GeV data sets, which are very 
similar. Our calculations are performed at Pbeam = 
286 GeV, since most of the data was obtained at this 
energy. We use the PDFs of the CBT model QUI at a 
fixed Q 2 of 25 GeV 2 , which is approximately the mean Q 2 
of the NA10 experiment. The top right panel shows the 
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FIG. 3: The ratio of Drell-Yan cross sections for " nv , , ."; 

0"-^ r (7T — +Au) 

(upper right) as a function of 



(upper left) and "w^-"!^"' 
X2, using nuclear PDFs from the CBT model for gold (red 
solid) and for N = Z nuclear matter (blue dot-dashed) . These 
calculations are performed for a pion beam energy of 160 GeV, 
at Q 2 = 25 GeV 2 . The lower panels show the same ratios but 
as a function of the dimuon mass, M, at a fixed Xi — 0.5. 
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ratio " Dv\ n -+m from the NA3 experiment [2 



The 

data was collected using a 150 GeV ir~ beam and the 
Q 2 range covered was 16.8 < Q 2 < 70.6 GeV 2 . Our 
calculations arc performed for Pbeam = 150 GeV and 
Q 2 = 25GeV 2 . 

The solid curves in Fig. 2 are calculations using the 
flavor-dependent Ua(%) and (Ia(x) from the CBT model 
with N/Z = 1.5, corresponding approximately to the 
N/Z values for the Au, W and Pt nuclei. The dashed 
curves correspond to the calculated ratios using the nu- 
clear PDFs from the CBT model for N = Z. Since 
ua/ud = dAJdo in this case, the dashed curves are 
representative of the predictions for flavor-independent 
EMC models. Figure 2 shows that the NA10 data do not 
exhibit a clear preference for the flavor-dependent ver- 
sus flavor-independent nuclear PDFs. In contrast, the 
NA3 data strongly favor the calculations using flavor- 
dependent nuclear PDFs. 

The Drell-Yan ratio, R±, of Eq. @ is an ideal exper- 
imental observable to search for flavor-dependent EMC 
effect, since it is directly proportional to cIa(x)/ua(x). 
The onl y ex isting data on R± is from the Omega collab- 
oration [29| ■ The lower panels of Fig. 2 show the ratios 



of 



"(tt^+W) 



collected with 39.5 GeV pions over a Q 2 



a DY (ir-+W) 

range of 5.3 < Q 2 < 7.3 GeV 2 and 16 < Q 2 < 25 GeV 2 . 
The calculations were performed at Q 2 = 7 GeV 2 and 



Q 2 = 20GeV 2 , respectively. Unfortunately, the existing 
data lack the precision for placing a useful constraint on 
the flavor dependence of the EMC effect. 

The plan for the COMPASS collaboration J3(| to mea- 
sure Drell-Yan cross sections with 160 GeV pion beams 
offers an exciting opportunity for a precise test of flavor- 
dependent EMC effect. Figure 3 shows the predictions 
of R± and R^ /D for the COMPASS kinematic coverage 
using the PDFs from the CBT model. The upper pan- 
els show the predictions as a function of x^ for a fixed 
dimuon mass of 5 GcV, while the calculations for a fixed 
x\ = 0.5 arc shown in the lower panels as a function 
of dimuon mass. The significant difference between the 
predicted ratios using the flavor-dependent versus flavor- 
independent nuclear PDFs provide a strong motivation 
for such measurements in the future. 

It is worth noting that several parametrizations of the 
nuclear parton distributions have been obtained from 
global analyses of DIS and Drell-Yan data |31h34{ . How- 
ever, pion-induccd Drell-Yan data were not included in 



these global analyses. The sensitivities of the pion- 
induced Drell-Yan data to the flavor dependence of the 
nuclear PDFs, as shown in this study, strongly suggest 
that these data should be included in future global anal- 
yses. 

In conclusion, we have examined the possibility of us- 
ing pion-induccd Drell-Yan processes as a sensitive exper- 
imental tool to study the flavor dependence of the modifi- 
cation of quark distributions in the nuclear medium. We 
suggest several Drell-Yan cross-section ratios as the most 
sensitive experimental quantities for such a study. The 
existing Drell-Yan data are not sufficiently accurate yet, 
although the NA3 data clearly favor the flavor-dependent 
over the flavor-independent nuclear PDFs. Precise future 
pion-induced Drell-Yan experiments can provide unique 
constraints that will help distinguish the various theoret- 
ical models and most importantly shed new light on the 
origins of the EMC effect. 

This work was supported in part by the U.S. Depart- 
ment of Energy and the National Science Foundation. 



[13 
[14 
[15 
[16 

[17 
[18 



[1] J. J. Aubert et al. (European Muon Collaboration), Phys. 

Lett. B 123, 275 (1983). 
[2] J. Gomez et al, Phys. Rev. D 49, 4348 (1994). 
[3] G. Bari et al., Phys. Lett. B 163, 282 (1985); A. C. Ben- 

venuti et al., Phys. Lett. B 189, 483 (1987). 
[4] A. M. Cooper et al, Phys. Lett. B 141, 133 (1984); 

H. Abramowicz et al., Z. Phys. C 25, 29 (1984). 
[5] K. Ackerstaff et al. (H ERMES Collaboration ), Phys. 

Lett. B 475, 386 (2000) [arXiv:hep-ex/9910071| . 
[6] J. Seely et al, Phys. Rev. Lett. 103, 202301 (2009) 

[arXiv:nucl-ex/0904.4448] . 
[7] P. Bordalo et al., Phys. Lett. B 193, 368 (1987). 
[8] D. M. Aide et al, Phys. Rev. Lett. 64, 2479 (1990). 
[9] P. L. McGaughey, J. M. Moss, and J. C. Peng, Annu. 

Rev. Nucl. Part. Sci. 49, 217 (1999). 
[10] D. F. Geesaman, K. Saito and A. W. Thomas, Ann. Rev. 

Nucl. Part. Sci. 45, 337 (1995). 
[11] P. R. Norton, Rept. Prog. Phys. 66, 1253 (2003). 
[12] M. Ericson, A. W. Thomas, Phys. Lett. B148, 191 

(1984). 

R. P. Bickerstaff, M. C. Birse, G A. Miller, Phys. Rev. 

Lett. 53, 2532-2535 (1984). 

I. G Cloet, W. Bentz and A. W. Thoma s, Phys. Lett. B 

642, 210 (2006) [arXiv:nucl-t h/060506l] . 

I. G Cloet, W. Bentz and A. W. Thomas, Phys. Rev. 

Lett. 102, 252301 (2009) [arXiv:nucl-th/0901.3559]. 

W. Bentz, I. G Cloet, J. T. Londergan and 



A. W. Thomas, arXiv:0908.3198 [nucl-th]. 

G P. Zeller et al. (NuTeV Collaboration) Phys. Rev. 

Lett. 88, 091802 (2002). 

I. Sick and D. Day, Phys. Lett. B 274, 16 (1992). 



[19] F. E. Close, R. L. Jaffe, R. G Roberts, and G G Ross, 

Phys. Rev. D 31, 1004 (1985). 
[20] G H. Llewellyn-Smith, Phys. Lett. B 128, 107 (1983). 
[21] M. Ericson and A. W. Thomas, Phys. Lett. B 128, 112 

(1983). 
[22] E. L. Berger, F. Coester, and R. B. Wiringa, Phys. Rev. 

D 29, 398 (1984). 
[23] I. G Cloet, W. Bentz and A. W. Thomas, forthcoming 

publication. 
[24] P. Souder et al, PR12-10-007, JLab Proposal to PAC-35, 

unpublished (2010). 
[25] B. Lu and B.-Q. Ma, Phys. Rev. C 74, 055202 (2006). 
[26] D. Dutta, D. Gaskell, K. Hafidi et al, PR12-09-004, JLab 

Proposal to PAC-34, unpublished (2009). 
[27] M. Gliick, E. Reya and A. Vogt, Z. Phys. C 53, 651 

(1992). 
[28] A. Michelini, CERN-EP/81-128 (1981); J. Badier et al., 

Phys. Lett. B 104, 335 (1981). 
[29] M. Corden et al, Phys. Lett. B 96, 417 (1980). 
[30] COMPASS Collaboratio, CERN Report No. CERN/SP 

SLC 96-14, SPSC/P 297. 
[31] K. J. Eskola, V. J. Kolhinen, and G A. Salgado, Eur. 

Phys. J. C 9, 61 (1999). 
[32] K. J. Eskola, H. Paukkunen 

04, 065 (2009). 
[33] M. Hirai, S. Kumano, and T. 

044905 (2004). 
[34] M. Hirai, S. Kumano, and T.-H. Nagai, Phys. Rev. C 76, 

065207 (2007). 



and G A. Salgado, JHEP 
H. Nagai, Phys. Rev. C 70, 



